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^^ I ' Abstract 

By going beyond Hubbard Hamiltonian we reflected correlation efi'ects 
accurately in the wavefunctions of H2- Using ab initio e-e interaction 
parameters resulted maximally entangled ground and third excited states. 
We assigned this maximally entangled character to the nonmagnetic (S=0) 
property of these states and also the minimally entangled character of the 

rf^ ■ first excited states to its magnetic property (87^ 0). By switching on a 

magnetic field an entangled state with Sz ~ can be extracted from a 
minimally entangled degenerate magnetic state. We suggest that presence 

J.— ^ ' of a moderate correlated system and a non-magnetic {Sz — 0) electronic 

state can be two criteria for finding maximally entangled electronic states 
in a realistic molecular system. 
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1 Introduction 

Since Einstein, Podolsky and Rosen [I] and Schrodinger [2] investigated the 
non-classical properties of quantum systems and entered new concept as en- 
tanglement in quantum physics, it had become strange property in interaction 
between particles. Recently the study of the entanglement is a useful resource 
for quantum communications and information processing [5] such as quantum 
teleportation [T, "F , superdense coding |5| , quantum key distribution [7^ , and 
quantum cryptography ^ whose input states are constructed to be maximally 
entangled. Also entanglement has been suggested as a quantitative measure for 
electron-electron (e-e) correlation in many body systems [HI HH] • As a simple 
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illustration of entanglement we can say that, if there is no way to write the 
states of two particles as a product of the two systems states in the Hilbert 
space, then there will be an entangled system [TT]. A lot of investigation has 
been done about measuring entanglement in the fermionic systems, such as the 
Woofers' measure [H] and the Schliemann's measure [131 HI]- Through Git- 
tings' investigation [TS], it is shown that all these entanglement measures are 
not suitable but the Zanardi's measure [16l [17] satisfies all desirable properties 
of entanglement measurement. 

The H2 molecule is the simplest two electron systems that can be used to imple- 
ment a robust many body calculation based on Hubbard model [IB] . Traditional 
Hubbard model which is a priory many-body approach usually is used as a first 
attempt to calculate entanglement. This model gives maximum entanglement 
by setting e-e interaction parameter [/ = 0, which is a controversial conclusion 
[rn [TOl HO] . In this paper we go beyond Hubbard model to calculate entangle- 
ment of the non-magnetic ground and magnetic excited states of H2 molecule. 
The H2 molecule is the simplest realistic two electron correlated system in nature 
which our model can be implemented and considering all direct and exchange 
interaction terms beyond Hubbard model let us account correlation effects more 
accurately in the wavefunctions. The Zanardi measurement for calculating en- 
tanglement was employed. We will investigate the ground and excited state of 
H2 to find the effect of the spin on the entanglement of states and results give 
maximally entangled states with non-zero U value. We also discuss the differ- 
ence between maximally, moderately and zero entangled states based on their 
magnetic property and correlation dependence. 

2 Calculation method 

Now we explain our calculation method. The complete Hubbard Hamiltonian 
is defined as [21] : 

^ = ^/^ iijCjg.CjCT + ^ / ^ ''ijhnC.i„Cjcr''^ma'Cla- \J-) 

ij(7 ijlraaa' 

The first term contains non- interacting part of Hamiltonian which can be written 
as: 

_^%Cj^Cjcr = £0 2_^ Cj^CjCT + t 2^ Cj^CjCT (2) 

ijc z— l,2;cr *7^j;cr 

^21 = tl2 ~ t 

where So is the energy of atomic orbital, cj^ and Cia are fermionic creation and 
annihilation operators on site i with spin a, respectively and t stands for the 
hopping integral between two H atomic sites of the electrons with the same a. 
The second term of Hamiltonian that contains e-e interaction part can also be 
written as: 
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where y^im contains all coulomb interaction between electrons which involves U 
as on-site Coulomb repulsion, / as inter-site Coulomb repulsion and nio- = Ci^Cia 
is density operator. The last terms Xi and X2 are the exchange interactions 
parameters that can be interpreted by quantum mechanics. We consider two 
electrons in two sites or Is orbital of H2 molecule with spin up and down, 
therefore we have four states for a single electron so for two electrons there are 
C(4, 2) = 6 states which are represented with notation ji/j >= \ni^ni^n2-\n2i > 
as 



m > 

\(P4 > 



11100 >,|(/52 > = 
10110 >,|(^5> = 



11010 >,|</73> = 
10101 >,|</76> = 



|1001 > 
10011 > 



(4) 



With these sets of states, we can write Hamiltonian as: 
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All the diagonal elements contain a term ta^ where it is roughly four times of 
the energy of an electron in the Is state of atomic hydrogen. Using ab initio 
energies for H2 and H2 from Ref. [22] and with the parametric energies of 
Hamiltonian, parameters £0, i, U^ J, Xi and X2, were evaluated and are given 
in Table I. 



Table 1: The calculated values of Hamiltonian parameters. 


Hamiltonian parameter 


£0 


t 


U 


J 


Xi 


X2 


calculated value (cV) 


-28.56 


6.36 


19.68 


17.90 


0.95 


1.36 



Entanglement measurement is defined by von Neumann's entropy as|16): 

S{pa) = -tr{pA log2 pa) (5) 

where A and B are bipartite subsystem which in our model are Is orbital of 
each Hydrogen atoms in our model. pA is reduced density matrix that is defined 
by: 

PA = trBPo = X! < ■?ls(|V' >< V'DIj >b 



Table 2: The parametric energies of hydrogen molecule are listed in the first 
column of the Table, where a±{x) = a; ± \/l + x^ and x = A ~„\ . The values 
of energies from Ref. [11\ are listed in the second column, eigenfunction for 
these states are also listed in the third column. The S'^ and 5*2 of each state are 
given in the next columns and the last column is the calculated entanglement 
S{pa) for these states, based on Zanardi's measurement. 
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where Itb stands for tracing over all sites except the B sites and | j >b is 
eigenstate for B part {\n2'\-n2i >) i.e. | 00 >, | 01 >, | 10 > and \ 11 >. 
After these calculation the reduced density matrix for the ground state (not 
normalized) becomes: 
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For other states, reduced density matrices have been evaluated accordingly and 
their entanglements were calculated. We listed the resultant entanglement val- 
ues in Table II. By acting S'^ on the \Eq),..., and ji^a) eigenstates, we could 
find total spin of each state in which \Eo), {£2) and [£^3) have S^ = 0, but \Ei) 
has S*^ = 2. We find eigenvalues and eigenstates (not normalized) of hydrogen 
molecule as shown in Table II, where a±{x) = x zL y/l + x^ and x = A^ ~^\ . 
Also the related Sz values of each eigenstates have been summarized in Table 
II. The eigenstates results schematically demonstrated in Fig. 1. Now let us 
discuss calculation results summarized in Table II and Fig.l. 



3 Results and discussion 

From this table we conclude that the many electron wavefunctions have very 
weaker dependence on the interaction parameters than the energy levels. The 
wavefunctions of the first and second excited states are explicitly independent 
from Hamiltonian parameters. The wavefunction of the ground and third ex- 
cited state are dependent upon e—e correlation parameters via a± = x±\/l + x^ 



with X 
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in Table II where the exchange interaction X2 is absent in the 



4(t+Xi) 

wavefunctions. Both of these states arc nonmagnetic (S=0). Using parameters 
of Table I and von Neumann's entropy of Eq. (6) an entanglement value S[p] ^2 
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Figure 1: The ground and excited states for H2 molecule before and after switch- 
ing on a magnetic field with their related eigenfunctions and spins. 



was obtained for these states. The maximum available value of entanglement is 
log2d for a system with the Hilbert space dimension of the smaller subsystem as 
d |23| . Accordingly, for H2 molecule the maximum available entanglement is 2 
therefore the resultant maximally entangled ground and third excited states in 
Table II can be explained by the corresponding wave functions of these states. 
The wavefunctions of the ground and third excited states are superposition of 
four body basis of the systems ipi^ (^3, ip^ and (pe with equal coefficients 
since the value of a± from table I becomes 1. 

From Table II, the first excited state Ei is a spin triplet state with 5 = 1 and its 
wavefunction is independent fromcorrelation parameters. The value of entangle- 
ment for the Sz = eignfunction is 1 where its value for the eignfunctions with 
iSz = ±1 is zero. The difference between entanglements of the degenerate wave- 
functions with different Sz can be explained by their related wavefunctions. The 
entanglement of the eignfunction oi \Ei, Sz = 0) state is a linear combination of 
ips and ip4 whereas the wave functions of other \Ei,Sz = ±1) are separable (1^92 
or (/J5). The importance of the nonzero spin of the first excited state is that this 
state can be detected by Electron paramagnetic resonance (EPR) under some 
condition [2j. When the magnetic field is off the wavefunction is a superpo- 
sition of the degenerate wavefunctions with different Sz values. We calculated 
entanglement of this degenerate wavefunction as 0. However, after switching on 
the magnetic field, a wavefunction with definite value of entanglement emerges. 
The values of the Sz and entanglement for this state are and 1 according to 
Table II. Therefore we conclude that in practice by switching on a magnetic 
field on a magnetic state, one can switch from a degenerate and not entangled 
wavefunction to an entangled wavefunction with Sz — Q- 

The second excited state E2 is nonmagnetic, 5 = 0, where the wavefunction 
is independent from correlation parameters. The calculated entanglement is 1. 
The wavefunction of this state is only a linear combination of two basis set Lp^ 
and (/3i, hence its entanglement is smaller than ground and third excited states. 
By comparing the values of the entanglement listed in Table II, we conclude 
that the nonmagnetic states (5 = 0) in which the wavefunction depend on in- 
teraction parameters are maximally entangled and the magnetic states whose 
wavefunctions are independent from correlation parameters are not entangled. 
Using von neumann's entropy in Eq. 6, we calculated the variation of entan- 
glement for maximally entangled ground state of H2 molecule with respect to 
the combination of correlation parameters x = ht+x\ ■ R-^sults were plotted in 
Fig. 2. 

In this figure we observe that the value of entanglement is maximum with 
X — Q {U — J). This conclusion resolves previously reported results from 
other groups [T71 [HI |20] who obtained maximally entangled ground state with 
U = Q. The [/ = result imply that the maximally entangled ground state is 
not attainable for the H2 since the ah initio value of U from Table I is 19.68 
which is very far from zero. In our given model the maximally entangled ground 
state obtained with U=J. This has a meaningful physical interpretation which 
states that in such physical system where the inter-atomic distance is very small, 




(U-J)/4(t+X ) 



Figure 2: Entanglement for the ground state of Hydrogen molecule versus e — e 
interaction parameters. 



the on-site Columbic repulsion U can be very close to the inter-site Columbic 
repulsion J. Indeed the H2 molecule is the best example of such systems when 
the inter-atomic distance is minimum or d = 0.7 A and the ab initio value of U 
and V from Table II gives x — 0.06. Using ab initio parameters of Table I, we 
obtained x = 0.06 and this value of x gives maximally entangled ground state. 
From Fig. 2 the maximum entanglement value for x = 0.06 is 1.99. 
As we observe in Fig. 2, in the extreme limit of | yj ^ 1, i.e. strongly correlated 
systems the entanglement becomes smaller and tends to 1. The non- magnetic 
property of this state, sets a limit on minimum available entanglement for this 
state. This point can be explained by using eigenfunction listed in Table II. 
For |a;| » 1 one of the a+ or a_ goes to zero and the other one becomes very 
large. In both cases the ground state wavefunctions listed in Table II reduce from 
extension on the four components to an extension on the two components similar 
to I £^2) state. Therefore the corresponding value of entanglement reduces from 2 
to 1. This also can be explained by tendency of the strongly correlated systems 
with I Y I ^ 1 to unpaired electronic configuration in atomic orbital such as (p2 
and ip5 states in which energy reduces by loosing U term in Hamiltonian. Such 
states tend to have parallel spin and magnetism with a reduced entanglement. 
In the case of if 2 where t, U and J are at the same order of magnitude (See Table 
I) and X is much close to zero hence the molecule is in moderately correlated 
regime and one can obtain the maximum available entanglement as it is shown in 
Fig. 2. In this regime we have both spatial and spin correlated wave function [17]. 
Neglecting exchange interaction in our model {Xi = X2 = 0), significantly 
displaces energy levels of the system (see Table II). However, the dependence of 
the ground and third excited state wavefunctions on the exchange parameters 
only is limited to Xi and putting Xi = yields x = 0.07 and similarly the 
maximum entanglement becomes 1.99 which is the same as the case of nonzero 
Xi {x = 0.06). In otherwords in order to obtain maximally entangled states, 
the most effective parameters are direct columbic interaction parameters U and 
J and exchange interaction parameters do not alter the value of entanglement 
significantly. 



4 conclusion 

In conclusion, in this paper we applied a robust many electron calculation on 
a simplest realistic two electron system i. e. H2 molecule. Going beyond tra- 
ditional Hubbard model let us to account correlation effects accurately in the 
many electron wavefunction of the ground and excited states. Using ab initio 
e-e interaction parameters, indicates to a moderately correlated regime for the 
molecule and resulted a maximally entangled ground and third excited state. 
The wavefunctions of the not magnetic (S=0) ground and excited states explic- 
itly depend on correlation parameters whereas the first excited states which is 
magnetic {S'^ — 2 and Sz 7^ 0) is not entangled. The second excited state is not 
magnetic but its wavefunction does not depend on correlation parameters there- 
fore it is a moderately entangled state. Anycase, by switching on a magnetic 



field an entangled state with Sz = can be extracted from a not entangled de- 
generate magnetic state. We suggest that in a realistic molecular scale systems, 
there is two criteria for finding maximally entangled electronic states, first the 
system should be in moderately correlated regime and second the system should 
have a non-magnetic {Sz = 0) electronic state. 
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